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Abstract 

Internal Diffusion Limited Aggregation (IDLA) is a model that describes the growth 
of a random aggregate of particles from the inside out. Shellef proved that IDLA 
processes on supercritical percolation clusters of integer-lattices fill Euclidean balls, 
with high probability. In this article, we complete the picture and prove a limit-shape 
theorem for IDLA on such percolation clusters, by providing the corresponding upper 
bound. 

The technique to prove upper bounds is new and robust: it only requires the ex- 
istence of a "good" lower bound. Specifically, this way of proving upper bounds on 
IDLA clusters is more suitable for random environments than previous ways, since it 
does not harness harmonic measure estimates. 

1 Introduction 

The Internal Diffusion Limited Aggregation (IDLA) model was introduced by Diaconis and 
Fulton in [7], and gives a protocol for recursively building a random aggregate of particles. 
At each step, the first vertex visited outside the current aggregate by a random walk 
started at the origin is added to the aggregate. This simulates the growth of an aggregate 
of particles from the inside out. 

In a number of settings, this model is known to have a deterministic limit-shape, mean- 
ing that a random aggregate with a large number of particles has a typical shape. On Z d , 
Lawler, Bramson and Griffeath [T3] were the first to identify this limit-shape, in the case 
of simple random walks, as an Euclidean ball. Their result was later sharpened by Lawler 
|12j . and was recently drastically improved with the simultaneous works of Asselah and 
Gaudillere [21 [3] and Jerison, Levine and Sheffield [101 [IDE]- O n other graphs, current 
knowledge is less precise. On groups with polynomial growth, the existence of the limiting 
shape is unknown, although Blachere gave bounds on the cluster [5]. On finitely generated 
groups with exponential growth, with a suitable metric, a limit-shape result was proved 
by Blachere and Brofferio [6]. Huss also studied IDLA for a large class of random walks on 
general non-amenable graphs in [8j. 

Another interesting question about IDLA is whether the limit-shape is robust to small 
perturbations of the underlying graph; For example, on the infinite cluster of supercritical 
percolation cluster of Z d . Shellef proved a sharp inner bound for the IDLA model on the 
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infinite cluster [13]. Figure [l] presents the IDLA aggregate built on the supercritical bond 
percolation cluster. 




Figure 1: IDLA aggregate with 6,300 particles on the supercritical percolation cluster on 
1? (edges are deleted w.p. 0.4)- Red points are in the aggregate, green points are in the 
cluster, and blue points are outside the cluster. 

To prove the existence of a limit-shape on the supercritical percolation cluster on Z d , 
it thus remains to show a sharp outer bound on IDLA. We provide such an outer bound, 
relying on Shellef's inner bound. 

Theorem 1. Fix d > andp > p c (7j d ) and let co be the infinite cluster of percolation on Z d 
with parameter p, conditioned on the origin o belonging to u. Let B {n) be the Euclidean 
ball of radius n centered at the origin. Let b (n) = \oj D B (n)\, the size of the part of the 
cluster that is in B (n). Let A bo ^(o) be IDLA generated by b Q (n) particles started at the 
origin. Then, for every e > and a.s. every oj, the IDLA on u) a.s. satisfies: 

B ((l — e)n) C Aw n )(o) C B {(1 + e)n) V n large enough. 

More generally, Theorem [2] below states that given a sharp inner bound on IDLA for a 
graph (with some regularity assumptions), a sharp outer bound exists as well. The inner 
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and outer bounds together prove the almost sure convergence to a limit-shape. The regular- 
ity assumptions are quite mild, and the theorem should be useful in great generality. Since 
the percolation cluster satisfies the necessary regularity assumptions, Theorem [2] implies 
Theorem [TJ We also mention other environments that satisfy the regularity assumptions, 
like Cayley graphs with polynomial growth (more details follow). 

Let us mention that, despite the fact that the upper bound on Z rf is not intrinsically 
harder than the lower bound, they invoke different ingredients. The lower bound usually 
harnesses estimates on the Green function, while the upper bound requires the use of 
upper bounds on harmonic measures. The Green function is fairly well understood in 
several (random) environments such as supercritical percolation, random conductances 
with elliptic condition, and so forth. This is not the case for harmonic measures. The 
techniques developed in this article allow to bypass this difficulty. Our main result relates 
an upper bound on IDLA directly to a lower bound. Roughly speaking, we show that if one 
knows that balls of a certain metric are contained in the aggregate, and not many particles 
are left over, then one can deduce an upper bound and shape theorem. 



1.1 Definition of IDLA 

Let G be a graph. Let S C G be a finite subset of G. In order to define IDLA, first define 
adding one particle started at x to an existing aggregate S. For x G G, denote by A(S; x) 
the IDLA aggregate obtained as follows: Let £ = (£(0),£(1), . . .) be a random walk on G 
started at £(0) = x and let ts be the first time this walk is not in S. Define 

A(S;x) :=SU{£(i 5 )}. 

It is standard to consider a slightly more general process, where the growth of the 
aggregate is stopped at certain stopping times, e.g., upon exiting a set T. Denote by 
A(S;x \— > T) the aggregate obtained by letting a particle randomly walk from x, but 
pausing it if it exits T, which is defined as follows. Let £ be a random walk on G started 
at x. Let ts be the first time this walk is not in S as above, and let tj< be the first time £ 
exits T. Define 

A{S;x^T) := S U {£(t s A (t T - 1))} . 
To keep track of the position of a paused particle, define 

P(S-x^T) = { £(* T-1 ) i£ Pr<ts, 
\ . otherwise, 

so _L means that the particle is already "absorbed" in the aggregate. 

Given vertices x\, . . . , x^ in G and a set T, define A(S; x\, . . . \-t T) to be the IDLA 
aggregate formed from an existing aggregate S by k particles, started at x\,. . . ,xt, and 
paused upon exiting T. That is, define inductively: So = S, Sj = A(Sj-\,Xj h-> T) for 
j £ {!,..., k}, and A(S; x\, . . . , x^ i->- T) = S^. Again, to keep track of paused particles, 
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define P(S; xi, . . . , x& i->- T) to be the sequence of particles paused in this process; Formally, 
if pj = P(Sj-i;xj i — y T) for j G {1, ...,&}, then P(5; xi, . . . ,Xfc h-> T) is the sequence 
fe' : When particles are not stopped, we define the aggregate similarly and we 

denote it by A(S; xi, . . . , x&). 

One reason to keep track of these paused particles is the so called Abelian property of 
IDLA: 

A{S; xi, ... , Xk) has the same distribution as 

A(A(S; xi,..., x k ^ T);P{S; xi, . . . , x k ^ T)). (1) 

Equation says that in order to sample A(S, x\,..., Xk), one can sample A(S; x\, . . . , Xk i-> 
T) while pausing particles upon exiting T and keeping track of them via P(S; x\, . . . , x^ h-> 
T), and then restart the paused particles on the obtained aggregate A(S; x%, . . . , X}~ >->■ T). 
For more details, see [71 [T3]. 

We are mostly interested in n particles starting at just one point. For an integer n > 0, 
by A n (x) we denote the IDLA aggregate with n particles started at x, that is, 

A n (x) = A^-x^^x). 

n times 

We also focus on pausing particles according to a metric p. Define the ball of radius 
r around x to be B x (r) = {y : p(x,y) < r}, and denote its size by b x {r) = \B x {r)\. As 
above, set 

A n (x i->- r) = A(0; x, . . . , x h-> S x (r)) 

and 

P n (x ^ r) = P(0; i 14 B x (r)). 

1.2 Assumptions 

We now make a few assumptions. The first two assumptions are independent of the IDLA 
process. We will always compare the IDLA cluster to a ball for a certain metric p on G, on 
which we now make a few hypotheses. Assume that p satisfies the following assumptions: 

Continuity (C): the metric is dominated by the graph distance do on G: there exists 
c > such that 

p(x,y) <c- d G (x,y) Vx,yGG. 
Regular volume growth (VG): there exist c, d > such that for every n > 0, 

< |S x (r)| < cr d Vx£ B (n) and n 1/d3 < r < n. 

The conditions on r and n allow to consider inhomogeneous environments. Keep in mind 
that the percolation cluster contains arbitrary finite graphs infinitely often. 
We will also make the following assumption on IDLA: 
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Weaker lower bound (wLB): there exists a > such that for every n > 0, 

F[B x (r) C A^/^x h-> r)] > a V x G S G (n + r) and n 1/d3 < r < n. 

In words, with noticeable probability, when releasing order b x {r) particles at x the aggregate 
contains B x (r). This assumption is easy to verify for many IDLA processes. In order to 
prove the existence of a limiting shape, we will assume a stronger lower bound for the 
aggregate grown around the origin, as we will see in the next section. 

1.3 Main results 

The first theorem relates a lower bound to an upper bound in well-behaved environments. 

Theorem 2. Let G be a graph and p be a metric on G satisfying conditions (C), (VG) 
and (wLB). Then, there exists a constant c\ > such that for any e > 0, 

r[A bo{n) (o) <£ B ((l + e)n) i.o] < P[| A bo(n) (o ^ n)\ < b (n)(l - Cl e d ) i.o.]. 

The previous theorem is especially useful when a lower bound is known. Indeed, al- 
though we need a statement stronger than a simple lower bound, the usual proofs of lower 
bounds always yield the fact that the IDLA process (almost) fills all large enough balls 
when stopped on their boundary: 

Lower bound (LB): \A bo r n \(o i— > n)\/b (n) converges almost surely to 1. 

The difference between (LB) and the "usual" notion of a lower bound is that the "usual" 
lower bound says that the aggregate eventually contains a large ball, whereas (LB) says 
that the aggregate contain a large ball, even when stopped upon exiting a (slightly bigger) 
ball. 

Corollary 3. Let G be a graph and p a metric on G satisfying conditions (C), (VG), 
(wLB), and (LB), then for every e > 0, a.s. 

J 4fe (n)(°) C B Q ((1 + e)n) V n large enough. 

Even though the main application of the theorem above and its corollary will be in the 
case of supercritical percolation, they apply in great generality and we believe that they 
will be useful in future works on IDLA. 

2 An upper bound on IDLA on percolation clusters 

Let us prove that Theorem [2] and its corollary imply Theorem [TJ In this case, p will be 
the Euclidean distance. By Theorem [2j we only need to check that (C), (VG), (wLB) and 
(LB) are satisfied: 
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(C) Since p <d w deterministically, property (C) is satisfied for every u with the constant 
1. 

(VG) Barlow proved in [3] that (VG) is satisfied for the distance d w for almost every 
environment. A classical result of pQ easily implies that for almost every environment 
u, there exists c\ = c\(uj) > such that for every n > 0, 

p(x,y) < d UJ (x,y) < c 1 p(x,y) V x, y £ B (n) : c\ log n < p(x, y) < n (2) 

so that (VG) is also satisfied for p with a possibly different constant (the result also 
follows from [1]). 

(wLB) In (T3|, Shellef proved that for any e > 0, there exists r] > such that the following 
holds for almost every environment cj: there exists C2 = 02(00) > such that 



[B ((l - e)n) C A bo[n) (o t-> n)] > 1 
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_C2_ 
d+2 



and P p [c2 > A] < e xv for all A > 0. All together, this implies that for almost every 
environment u, there exists C3 = 03(0;) > such that for every n > 0, 



B x ((l - e)r) C A bx{r) (x i-> r) VsG £ (n) and n 1 ^ 3 < r < n > 1 - — . (3) 

The condition (wLB) follows readily. The result in Shellef deals with the event 
B Q ((1 — e)n) C A bo ( n }(o) (particles are not stopped at distance n), but the proof 
actually implies this stronger result. 

(LB) Finally, note that the comparison between distances Q and (VG) imply that 

b (n) - b a ((l - e)n) < c 4 e d b (n) 

for some constant C4 = c^uj) > depending on the environment, so that ^ implies 
(LB) for almost every environment. 



3 Proof of Theorem [2] 

From now on, we fix a graph G satisfying (C), (VG) and (wLB). Constants in the proof 
always depend only on the constants involved in (C), (VG) and (wLB), i.e. c,d and a. 

The following lemma shows that a lower bound on the aggregate implies a lower bound 
on hitting probabilities. It is a general statement not invoking any of the conditions (C), 
(VG) or (wLB). In the following, we make a slight abuse of notations: £ will denote a 
random walk as well as its trace. 



6 



Lemma 4. Let Q C B C G and x £ B. Let £ be a random walk started at x and stopped 
on exiting B. For any t > 0, 

P[C n Q ^ 0] > P[B c A t (x ^ B)] ■ \Q\/t. 

The above lemma is most useful when t is chosen so that ¥[B C At(x \— > B)] is of order 1. 

Proof. Let £i, . . . , £t be the t independent random walks started at x and stopped on exiting 
B that generate the aggregate At(x \— > B). Let J G {1, . . . , t} be a uniformly chosen index 
independent of the random walks. Consider the set T of j G {1, . . . ,t} so that £j hits Q 
before exiting B. Since Q C B, the inclusion B C Ajfi 4 fi) implies |T| > \Q\. Since J is 
independent of £i, . . . , £t, 

p[^j n Q / | b c ^(x ^B)]>P[Jer|Bc A t (x t-> b)] > |Q|/t. 

The lemma follows since the distribution of £j is that of a random walk started at x and 
stopped when exiting B. □ 

By assumption on G, we can hence get the following hitting probability estimate, 
which states that a random walk hits a set, whose complement has size at most er d , with 
a probability that is bounded away from zero. 

Lemma 5. There exist e, r/ > such that for large enough n and n 1 ' < r < n, the 
following holds. Let x G B (n) and let S C B (n + r) be so that \S \ B Q (n)| < er d . Lei £ 
be a random walk started at x and stopped upon exiting B a (n + r). Then, 

F [C n (B n (n + r)\(SU B (n))) ^ 0] > 77. 

Proof. For every path 7 from inside B {n) to outside B Q {n + r), let 2/(7) be the first vertex 
on 7 so that p{y{"f),B (n)) > r/2. Denote by Y the set of all 2/(7) for such paths 7. Every 
path from x to outside B Q (n + r) must hit Y. By Markov's property, it thus suffices to 
prove the theorem for starting points y £Y. Fix y eY. 

Let B = B y {r/2>) and Q = B\S. By iVG) and by assumption on 5, 

|Q| > I( r/3 )rf_ er rf> i(r/4) d , 
c c 

with e = 4" d /c. By (wLB) with i = 6<,(r/(3a)), 

P[J3 C A t {y i-> r/3)] > a. 

Let £ be a random walk started at y and stopped on exiting B. Lemma [4] and (VG) imply 
that 
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Note that r/2 < p(y, B (n)) < r/2 + c thanks to the definition of Y and (C). Therefore, 
^(y, G \ B (n + r)) > r — r/2 — c > r/3 for n large, and so B C B Q (n + r) \ B (n). We 
deduce 

P [e n (B„(n + r)\(SU S (n))) / 0] > P[£ n Q / 0] > r,. 

□ 

After analyzing the behavior of a single particle, we can analyze the behavior of the 
whole aggregate. The following lemma says that, with high probability, a constant fraction 
of the aggregate is absorbed in a wide enough (yet still very fine) annulus. 

Lemma 6. There exist 5 > and p < 1 such that for all n large enough, for all n l ^ dJr1 ^ < 
k < n and xi, . . . ,Xk £ B (n), and for all S C B {n), 

¥[\A{S ]Xl ,...,x k ^ B {n + k 1 ' d ))\S\ <5k] <p k . 

Proof. Let r = k 1 ^. Fix e, rj as in Lemma[5] Let £i, . . . , ^ be the random walks started at 
x\, . . . , Xk that generates the aggregate. Let k' = [ek\ < er d . For j £ {1, . . . , k'}, denote 

Aj = A(S; xi, . . . , Xj i— > B a (n + r)). 

Since \Aj \ B (n)\ < j < er d , Lemma [5] implies that for all j £ {1, . . . , k'}, 

P fe+i n (B (n + r) \ Aj) ± | A 3 \ > V . 

Therefore, \A(S; x±, . . . , x^ h-> i? D (n + r)) \ 5| dominates a (A/, r/)-binomial random variable. 
Thus, there exist 5 > and p < 1 depending only on e, 7/ such that 

F[\A(S; Xl ,...,x k ^ B a (n + r)) \ 5| < <5fc] < p fc . 

□ 

We now turn to the proof of Theorem [2} The proof consists of inductively constructing 
a sequence of aggregates Aj by pausing the particles at different distances rij from the 
origin. If kj is the number of paused particles, we choose the next distance Tij+i, at which 
we pause the particles again, in terms of nj and kj. We iterate this procedure until there 
are less than n 1 '^ +1 ^ paused particles. At this point, there are too few particles to matter. 

Proof of Theorem^ Fix n > 0. Define Aj,nj,Pj, kj as follows: 

• Let no = n and Aq = A bo ^{o h-> B (n)). Let Pq = P bo ^(o \-t B (n)) and let 
k = \P \. 
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For j > 0, define 

' n 3 +k) /d if kj >n 1 /(^+D, 



n j+ i 



oo otherwise. 



Let Aj+i = A(Aj;Pj i->- B (nj+i)). Let Pj + \ = P(Aj-,Pj i-> B (nj+i)) and let 
= l-Pj+il- 

Let J be the (random) first time at which kj < n l K d+l \ By construction, Aj = Aj + \ for 
any j > J + 1. The Abelian property ([I]) guarantees that Aj + \ and .A& o ( n )(o) have the 
same law. 

By construction, Aj C B (nj). Since fcj < n 1 /^ 1 ) and p is continuous (C), the fcj 
last particles cannot grow long arms. Formally, Aj + \ C B (nj + cra 1 /^ 1 )). 
Since J < n, by Lemma[5j for some 5 = 5(a, c, d) < 1, 

F[3 1 < j < J : kj > (1 - 5) j k } < F[3 1 < j < J : kj > (1 - <5)%_i] < np nl/(d+1) . 
This implies that with probability at least 1 — np nl/td+1) , 

nj = n + k^ d H h k 1 /^ < n + ■ 



and if nj + Cn 1 /^ 1 ) > (1 + e)n, then > (en - Cn 1 /( d+1 ))(l - (1 - 5) l / d ). So, using 
(VG), for any e > 0, 

{A 6o( „)(o) £ B ((l + e)n)} C {^o > ci£ d 6 (n)} U {3 1 < j < J : kj > (1 - . 
Using the Borel-Cantelli Lemma, we deduce the result easily. □ 
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